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Abstract 

The report reviews the problem of topological coherent modes, which 
are nonlinear collective states of Bose-condensed atoms. Such modes 
can be generated by means of alternating external fields, whose frequen- 
cies are in resonance with the transition frequencies between the related 
modes. The Bose gas with generated topological coherent modes is a 
collective nonlinear analog of a resonant atom. Such systems exhibit a 
variety of nontrivial effects, e.g.: interference fringes, interference current, 
mode locking, dynamic transitions, critical phenomena, chaotic motion, 
harmonic generation, parametric conversion, atomic squeezing, and en- 
tanglement production. 



1 Introduction 

Trapped Bose gas, being dilute and cooled down to temperatures below the Bose- 
Einstein condensation temperature, forms a coherent system well described by the 
Gross-Pitaevskii equation [1-5] . Stationary solutions to the latter are called the topo- 
logical coherent modes [6]. The ground state, corresponding to the lowest energy level 
of the Gross-Pitaevskii stationary equation, represents the standard Bose-Einstein 
condensate. The states, associated with the higher energy levels, describe nonground- 
state condensates. The higher topological coherent modes can be generated by means 



of alternating fields, with the frequencies that are in resonance with the chosen tran- 
sition frequencies [6]. 

The name topological coherent modes comes from the following. Different stationary 
solutions to the Gross-Pitaevskii equation, related to distinct energy levels, possess 
essentially dissimilar spatial shapes, with different number of zeros. The modes, dis- 
playing distinct spatial topology, can be named topological. These modes can be 
characterized as coherent due to the fact that the Gross-Pitaevskii equation can be 
interpreted as an exact equation for coherent states [7]. 

The topological coherent modes, described by the nonlinear Gross-Pitaevskii equa- 
tion, should not be confused with elementary collective excitations, defined by the 
linear Bogolubov - de Gennes equations. The elementary excitations describe small 
oscillations around a given nonlinear topological mode and do not change the topology 
of the latter [1]. 

The general notion of the nonlinear coherent modes and the possibility of their 
generation by means of resonant fields was advanced in Ref. [6]. A particular case 
of vortex creation was considered in Refs. [8,9]. Various properties of topological 
coherent modes were also studied in theoretical papers [10-27]. A dipole topological 
mode was generated in experiment [28]. The feasibility of resonant formation of multi- 
mode condensates, consisting of several topological coherent modes, was investigated 
[29,30]. 

Bose gas of trapped atoms with resonantly generated topological coherent modes 
exhibits a variety of unusual features. The aim of this report is to present a general 
picture of such a resonant system (Section 2) and to give a survey of its most interesting 
properties (Section 3). 



2 Topological Coherent Modes 

The topological coherent modes are defined [6] as the stationary solutions to the 
eigenvalue problem 

H[(p n ] (p n (r) = E n <p n (r) , (1) 

which is represented by the stationary Gross-Pitaevskii equation with the nonlinear 
Hamiltonian 

HM = -^ 2 + U{v)+NAM\ (2) 

containing a trapping potential U{v) and the interaction intensity A s = A-K% 2 a s /m, 
with m being atomic mass; a s scattering length; and N the total number of atoms. 

In what follows, we shall use the notation for the scalar product (ipm,<Pn) = 
J (fi m (r)<p n (r)dr. The eigenfunctions (fi n {r), because of the nonlinearity of problem (1), 
are not necessarily orthogonal, so that (v?m, <fin) is not compulsory the Kroneker delta 
5mn- But the functions <Pn(r) can always be normalized, with ((p n ,(p n ) = ||^n|| 2 = 1. 

If the trapped Bose gas is initially in one of the coherent modes n, then to generate 
another mode requires applying an alternating field 

V(r, t) — Vi(r) cosujt + V2(r) sin cot , 

with a frequency ui being close to one of the transition frequencies 

Umn = -T {E m — E n ) . (3) 
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Say the modes with the energies Ei and E2 are connected, such that E\ < E2, with 
the transition frequency being u>2i- The related resonance condition reads as 

AW -.A / A\ 

< 1 , Au) = UJ — LU2\ ■ (4) 

With an applied time-dependent field, we have the temporal Gross-Pitaevskii equa- 
tion 

ih ^ <p(r,t) = (H[<p] + V) <p(r,t) , (5) 
in which V = V(r, t). The alternating field can be represented as 

V(r,t) = is(r)e- t + is*(r) e -- t , (6) 

where B(r) = Vi(r) — iV2(r). Experimentally, the alternating field can be realized as 
the modulation of the trapping potential. Another way could be by alternating the 
scattering length a s (t) by modulating an external magnetic field close to a Feshbach 
resonance [31-33]. 

One possibility of studying the resonant generation of coherent modes would be 
by a direct numerical solution of the temporal equation (5), which we have done 
in our works [29,30]. However to get a deep physical insight into the problem, it 
is necessary to develop an analytical theory. Such a general theory, based on the 
averaging technique [34], was developed in Refs. [6,12,15,20-24,29,30]. In order to 
show that the basic equations of the analytical approach can be accurately derived 
and all consideration is well mathematically grounded, we describe, first of all, the 
main steps of this derivation. 

We can look for the solution of Eq. (5) in the form of the mode expansion 

<p(r, t)=Y^ c n {t)>fi n {r) exp (- l - E n tj , (7) 

n 

with the coefficients c n (t) being slow functions of time, such that 

do, 



E n 



dt 



« 1 . (8) 



Then the functions c n (t) can be treated as temporal quasi-invariants with respect 
to the fast exponentials exp(—iE n t/h). Substituting expansion (7) into Eq. (5), we 
multiply the latter by the mentioned exponential and average over time according to 
the rule 



1 f T 

hm — / f(c„,t)dt 
T ^°° 1 Jo 



where the quasi-invariants c n are kept fixed. Averaging in this way the equality 
||v?|| 2 = 1, we have the normalization condition 

^>„(*)| 2 = 1. (9) 

n 

For what follows, we need the notation for the transition amplitudes: one of them 
being due to atomic interactions, 

a mn = A 3 ^ (|^m| 2 , 2|<^ n | 2 - |^ m | 2 ) , (10) 
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and another related to the external field (6) , 

[3mn = i (p m , B<fin) , (11) 

where B = B(r). Note that from normalization (9) one has jc„| 2 = 1 — ^ m (^n) \ Cm \ 2 ■ 
Using this, the nonlinear part, resulting from Eq. (5), can be represented as 

V — ^ 2 2 X — ^ 2 

A s — y [2 — Jrrm)(|ln| , |</>m| )= anm|c m | C n + «„„C„ . 

Thus we come to the equation 

. dc n \ ^ I 1 2 . 1 r ^ iAui-t 1 r ,o* — iAiu-t /in\ 

1 = 2-*i anm ' Cm ' Cn + Q n» c » + 2 " nl " 12C2e + 2 " n2 " 12Cie ■ y 12 > 

What is of physical interest is the behaviour of the fractional mode populations 

w n (t) = \c n (t)\ 2 . (13) 

These do not depend on the phase of c„. Therefore, we may employ the gauge 
transformation 

c„ — > c n cxp(-ia nn t) . (14) 

Then Eq. (12) reduces to 



* = a ™ m \ Cm \ 2cn + 2~ ^iPi2C2e' A ' t + i <5„2/3i2Cie lA '* , (15) 



where 

A = Auj + an — «22 • (16) 

For neighboring modes, one has an m 022 and A w Aw. Otherwise, it is always 
possible to choose such a detuning Au> that A be small, which is assumed in what 
follows, A I < US. 

To preserve well-defined resonance, it is necessary, in analogy with the case of 
resonant atoms [35], that the transition amplitudes, involved in the process, be small. 
In the considered problem, such amplitudes are given by Eqs. (10) and (11). So, it is 
necessary that the transition amplitudes, due to atomic interactions, be small, 



Q12 



«1, 



Ct21 



W21 



« 1 , (17) 



W21 

as well as the amplitude determined by the modulating field, 

^ « 1 . (18) 

U)2\ 

From Eq. (15) it stems that if at the initial time c„(0) = for n ^ 1, 2, then always 

c„(t) = 0, (n/1,2) (19) 
for all t > 0. Hence Eq. (15) reduces to the system of two complex- valued equations 

• 1 2 . 1 r, iA-t • dC2 I 1 2 . 1 r,* -iA-t ; on \ 

1 -^j- = ai2\C2\ ci + - pi2C 2 e , i — = c*2i|ci| c 2 + - /5i 2 cie . (20) 
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Such equations, though resembling the case of two coupled electromagnetic modes [36] , 
differ from that by the presence of the nonlinearity caused by atomic interactions. 

The complex-valued system (20) is equivalent to four real-valued equations. How- 
ever, because of the global gauge symmetry and due to the normalization condition 
ci | 2 + |c2 1 2 = 1, the related dynamical system is, actually, two-dimensional. To show 
this, it is convenient to define Cj = \cj\ exp(wr;t) and /3i2 = /3e 17 , where f3 = \(3i2\- Let 
us also introduce the notation 

a = 2 ( ai2 + a21 ) ' ^ = ^ + 2 ' ai2 _ a ' 21 ^ ' 
Then for the population difference 

s = M 2 -M 2 (2i) 

and effective phase difference 

aj = 7r 2 -7ri+7 + A, (22) 
we find the system of two equations 

ds _ r dx f3s 

— = — pv 1 — smx , — — as H ■ cos x + d . 23) 

dt dt VI - s 2 

The consideration can be generalized to the case of the multiple generation of 
topological coherent modes [29,30]. This requires, instead of one modulating field (6), 
the action of several oscillating fields 

V(r, t) = lYl [ B ^y^ + ^We-^'] , (24) 

3 

whose frequencies are tuned to the resonance with the chosen transition frequencies. 
For instance, in the case of three coexisting modes, we obtain 

i = (ai2|c 2 | 2 + ai3|c3| 2 ) ci + fi , 
i = (a 2 i |ci| 2 + a 2 3|c3| 2 ) c 2 + fc , 

i ^ = |ci | 2 + a 3 2|c 2 | 2 ) a + f A , (25) 

where the forces fj, related to the modulating fields, depend on the type of the 
mode generation, whether this is the cascade generation, V-type, or A-type generation 
[29,30]. The effective detuning A mn = Aw m „ + a„„ — a mm , where Alo 
is again assumed to be small, |A m7l | <JC uj mn . 

In the same manner, one can derive the evolution equations for the amplitudes c n (t) 
for an arbitrary number of generated coherent modes, whose populations are given 
by c n | . The resulting equations, such as (20), (23), or (25), are nonlinear because 
of the binary atomic interactions. One could also include three-body interactions, 
which would yield the fifth-order nonlinearity with respect to |c„|. Such three-body 
interactions can play an important role in describing the dissipation caused by three- 
body recombinations [1,37]. 

One could also take into account nonadiabatic corrections to the atomic evolution 
equations. However nonadiabatic description is of vital importance only for noncon- 
fined motion of atoms [38] , while for trapped atoms nonadiabatic corrections amount 
to at most a few percent [39]. 
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3 Dynamic Resonant Effects 



Trapped condensed Bose gas, with resonantly generated topological coherent modes, 
resembles a resonant atom, hence, such a resonant condensate should display the 
features typical of resonant finite-level atoms. But this resonant condensate is, in 
addition, a collective nonlinear system, because of which it can possess many other 
unusual properties, not existing in finite-level atoms. 

(1) Interference Fringes. The total density p(r, t) = \<p(r,t)\ of trapped atoms, 
with generated topological modes, is not simply the sum of the partial mode densities 
p n (r, t) — |c„(i)v?„(r)| 2 , but the interference fringes arise, described by the interference 
density 

/w(r, t) = p(r, t) - ^ p„(r, t) , 

n 

with fast oscillation in time [20,23]. 

(2) Interference Current. Similarly, there exists a fastly oscillating interference 
current 

jint(r,t) =j(r,t) - y^j n (r, t) , 

n 

sometimes called the internal Josephson current [20,23]. 

(3) Mode Locking. Under this effect, the fractional mode populations w n are 
locked in the vicinity of their initial values, so that either 

< w n (t) < i 

for all t > 0, or 

\<w n {t) < 1 , 

never crossing the line w n — 1/2, but being either below it or above it, depending on 
initial conditions [6,23,24]. 

(4) Dynamic Transition. Varying the system parameters, the dynamics of the 
mode populations can be qualitatively changed from the mode locked regime to the 
mode unlocked regime, when the mode populations fluctuate in the whole region 

< w n (t) < 1 , 

independently of their initial values [15,23,24]. 

(5) Critical Phenomena. On the manifold of the system parameters, there exists 
a critical surface in the vicinity of which the dynamics of the mode populations be- 
comes unstable. Crossing the critical surface, when varying some parameters, changes 
the dynamics between the mode locked and mode unlocked regimes. Close to the sur- 
face, a tiny variation of some of the system parameters, say of the pumping amplitude 
or of the detuning, provokes drastic changes in the dynamics of the mode populations 
[15,16,20,23,24]. The location of the critical surface also depends on the initial setup. 
Thus, for the case of two coexisting topological modes, we have the critical surface 
described by the relation 

i qsq — /3^/ 1 — Sq cosa^o + Ssq = (3 sgn a , 



G 



in which so = s(0) and xo = x(0). For each given initial conditions, this is a surface in 
the three-dimensional space of the parameters a,/3, and 8. For so = =Fl, the critical 
surface reads as 

(3 sgn a ± 8 = - a . 

Fixing so = — 1, we can reduce the above relation to two critical lines on the manifold 
of the parameters 

N ' H ' 

These lines are 

6 + (Q>0), 

A time-averaged system displays on the critical surface critical phenomena typical of 
statistical systems with the second-order phase transitions. For the averaged system 
it is possible to define an effective capacity and susceptibility, which diverge on the 
critical surface [15,16,20]. 

(6) Chaotic Motion. Fractional mode populations for a two-mode condensate are 
always periodic functions of time. But if the number of modes in the condensate is 
three or larger, then, depending on the system parameters, one has either quasiperiodic 
or chaotic motion. For instance, in the case of the three-mode condensate, for which 
otij = a, Pij — [3, and A;j = 0, chaotic motion appears when 



(a < 0) 



> 0.639448 , 



a 

that is, under a sufficiently strong pumping [29,30]. 

(7) Harmonic Generation. The generation of topological coherent modes may 
occur not solely under the direct resonance condition lo — u>2i, but also under the 
condition of harmonic generation 

ncu = UJ21 (n = 1, 2, . . .) , 

when just one modulating field, with a frequency ui, is applied [29,30]. 

(8) Parametric Conversion. Another possibility of generating the topological 
modes is under the condition of parametric conversion 



^(±Wj) = W21 , 



when several alternating fields, with frequencies u)j, are involved [29,30]. 

(9) Atomic Squeezing. Squeezing is a quantum effect that does not exist for 
classical quantities. So, for treating it, one has to quantize the modes by consid- 
ering the mode amplitudes c n as Bose operators. Then employing the pseudospin 
representation, one can define the squeezing factor 

n- 2A 2 (5.) 
Q ~ \< 5± > | ' 
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where A 2 (S Z ) is the dispersion of the z-component of the total spin and S± arc 
rising and lowering operators, respectively. The defined squeezing factor describes 
the relation between the dispersion, associated with the mode populations, and the 
dispersion of a relative phase. When Q < 1, one says that the atomic squeezing occurs. 
Then the mode populations can be measured with a higher accuracy as compared to 
the measurement of current [23]. 

(10) Entanglement Production. In the condensate with topological coherent 
modes, entanglement can be produced. A general measure of entanglement produc- 
tion, valid for arbitrary systems, was introduced in Ref. [40]. The measure of entan- 
glement, generated by an operator A, is defined as 

£ (^log 



\A®\\t> ' 

where || ■ ||x> implies the restricted norm over a set V of disentangled functions, and 
A® is a nonentangling counterpart of A. Considering a p-particle reduced density 
operator p p for a multimode condensate, we find 

e(p p ) = (1 - p) log sup w n , 

n 

where w n = w„(t) = |c„(t)| 2 . Hence the above measure e(p p ) = £ P (t) is a function of 
time, characterizing the entanglement evolution [40]. 



Summarizing, we have shown that by applying resonant modulating fields to a 
Bose-condensed gas of trapped atoms, one can create topological coherent modes. We 
have developed an analytical theory describing the condensate with such coherent 
modes and also made numerical calculations by directly solving the Gross-Pitaevskii 
equation [29,30]. All results of the analytical theory are in very good agreement 
with numerical solutions. The condensate with topological coherent modes is a novel 
system possessing a rich variety of nontrivial properties which could be employed in 
many applications. 
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